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ABSTRACT 

We observe a sharp transition from a singular, high-mass mode of star forma- 
tion, to a low-mass dominated mode, in numerical simulations, at a metallicity 
of 10~^ Zq. We incorporate a new method for including the radiative cooling 
from metals into adaptive mesh-refinement hydrodynamic simulations. Our re- 
sults illustrate how metals, produced by the first stars, led to a transition from 
the high-mass star formation mode of Pop III stars, to the low-mass mode that 
dominates today. We ran hydrodynamic simulations with cosmological initial 
conditions in the standard ACDM model, with metallicities, from zero to 10~^ 
Zq, beginnning at redshift, z = 99. The simulations were run until a dense core 
forms at the center of a 5 x 10^ Mq dark matter halo, at z ~ 18. Analysis of the 
central 1 core reveals that the two simulations with the lowest metallicities, 
Z = and 10~^ Z©, contain one clump with 99% of the mass, while the two with 
metallicities, Z = 10~^ and 10~^ Zq, each contain two clumps that share most of 
the mass. The Z = 10~^ Zq simulation also produced two low-mass proto-stellar 
objects with masses between lO'^ and 10"^ Mq. Gas with Z > lO'^ Zq is able 
to cool to the temperature of the CMB, which sets a lower limit to the minimum 
fragmentation mass. This suggests that the second generation stars produced a 
spectrum of lower mass stars, but were still more massive on average than stars 
formed in the local universe. 

Subject headings: stars: formation 



1. Introduction 



Numerical simulations have shown that the very first stars invariably formed in isolation 
and were much more massive than the sun, due mainly to the inability of primordial gas 
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to effi c iently cool at low temp eratures (lAbel et al.l 120021 : iBromm et al.l l2002l : lYoshida et al. 
20061 ). iTumlinson et al.l (120041 ) have suggested that the Pop III IMF was not dominated by 
very massive stars (M > 140 Mq), but instead by stars with M = 8-40 Mq. Even this IMF, 
though, is still remarkabl y distinct from that observed for the local universe , which peaks at 
less than one solar mass (IMiller fc Scald Il979l : iKroupal l2002l : IChabrieii 120031 ) . 



The deaths of the first stars produced and distributed copious amounts of metals into 
their surrounding s, through either core-c ollapse (M ^ 10 Mq) or pair-instability (M ^ 140 
Mq) supernovae (iHeger fc Woosleyll2002l ). These metals provide additional avenues for ra- 
diative cooling of the ambient gas, through fine-structure and molecular transitions, as well 
as continuum emission from dust formed from the supernova ejecta, permitting the gas that 
will form the next generation of stars to reach temperatures lower than what is possible for 



metal-free gas. Fragmentation of collapsing gas wi 
decreasing in temperature as the density increases ( 



1 continue so long as the gas can keep 



Larson 



20051 ) ■ or until the gas becomes 



optically thick to its own emission (ILow fc Lynden-Belllll976l ). The minimum fragment mass 
is determined by the local Jeans mass. 



Mj ~ 700 MQ{T/200Kf/^{n/10'^cm 



W2)- 



(1) 



where T, n, and /i are the teniperat ure, number density, and mean molecular weight, at the 
halt of fragmentation (jLarsonll2005l ). For metal- free gas, a minimum temperature of ~ 200 



K is reached at n ~ 10^ cm ^ when coolirig bec omes inefficient, yielding a Jeans mass 



Mj ~ 10^ Mq (lAbel et al.l 120021 : iBromm et al.ll2002l ). At some certain chemical abundance, 
it is conjectured that metals provide sufficient cooling, so that the temperature of the gas 
continues to decrease as the density increases past the stalling point for metal-free gas, 
allowing the collapsing gas-cloud to undergo fragmentation and form smaller and smaller 
clumps. The enrichment of gas to some critical metallicity, Z^r, will trigger the formation of 
the first low-mass (Pop II) stars in the universe, as the gas can cool to lower temperatures at 
higher metallicity, in general. The value of Zcr can be estimated by calculating the metallicity 
required to produce a cooling rate equal to the rate of adiabatic compression heating at a 



given temperat ure and density. 



his h as been carried out for individual alpha elements, such 



as C and 0. bvlBromm &: Loebl (120031 ) . and C, O, Si, Fe, as well as solar abundance patterns 



bv lSantoro fc Shulll J2006h . yielding roughly, lO-^-^ Zq ^ Z^r ^ lO^^ Z©. 



Aside from the minimum clump mass, however, not mu ch more can b e said about the 
spectrum of clump masses produced during fragmentation. lOmukai et al.l (120051 ) use one- 
zone models with very sophisticated chemical networks to follow the evolution of temperature 
and density in the center of a collapsing gas cloud, for a range of metallicities. The predic- 
tions of fragmentation from this work, though, are based solely on statistical arguments of 
elongation in prestellar cores and do not capture the complex processes of interaction and ac- 
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cretion associated with the formation of multiple stars (iBate et al.ll2003l ). iTsuribe fc Omukai 



( 120061 ) simulate the high density (n > 10 cm ) evolution of extremely low-metallicity gas 



(Z < 10~^ Zq), but the conclusions of this work are limited by the fact that the simulations 
are initialized at an e xtremely late phase in the evolution of the prestellar core. The nu- 
merical simulations by lBromm et al.l (120011 ). which use cosmological initial conditions, show 
fragmentation in gas with Z = 10~^ Z©, but a mass resolution of 100 Mq prevents this study 
from saying anything conclusive about the formation of sub-stellar mass objects. 

In this paper, we present the results of three-dimensional hydrodynamic sir nulations of 
metal-enriched star-formation. These simulations are similar in nature to those of Bromm et al 



(I2OOII ). but with vastly improved numerical methods and updated physics. We describe the 
setup of our simulations in §21 with the results in §3] and a discussion of the consequences of 
this work in ® 



2. Simulation Setup 



We perform a series of four simulations, with constant metallicities, Z = (metal-free), 
10~^ Zq, 10~^ Zq, and 10~^ Z^,, using the Eulerian adaptive me sh refinement hydrodynamics/N- 
body code, Enzo (IBryan fc Norman! 119971 : lO'Shea et al.ll2004l ). The metallicity is held con- 
stant throughout each simulation in order to isolate the role of heavy element concentration 
in altering the dynamics of collapse compared to the identical metal-free case. In reality, 
metals will be injected over time into star forming gas by Pop III supernova blast waves, 
and the mixing of those metals with the gas will not be completely uniform. Here we fo- 
cus on an idealized approximation in order to capture the essential physics of collapse and 
fragmentation. 

Each simulation begins at z = 99, in a cube, 300 h~^ kpc comoving per side, in a ACDM 
universe, with the following cosmological parameters: Qm = 0.3, Q\ = 0.7, Qb = 0.04, and 



Hubble constant, h = 0.7, in units of 100 km s ^ Mpc ^. We init ialize all the simu . 



ations 



identically, with a power spectrum of density fluctuations given by lEisenstein &: Hul (119991 ) . 
with as = 0.9 and n = 1. The computation box consists of a top grid, with 128^ cells, and 
three static subgrids, refining by a factor of 2 each. This gives the central refined region, 
which is 1/64 the total computational volume, an effective top grid resolution of 1024'^ cells. 
The grid is centered on the location of a ~ 5 x 10^ Mq dark matter halo that is observed 
to form a t z ^ 18 in a prior dark-matter-only simulation, as is done similarly in lAbel et al. 



( 120021 ): lO'Shea et al.l (120051 ). Refinement occurs during the simulations whenever the gas, 
or dark matter, density is greater than the mean density by a factor of 4, or 8, respectively. 
We also require that the local Jeans length be resolved by at least 16 grid cells at all times 
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in order to avoid artificial fragmentation as prescribed by iTruelove et al.l (119971 ) 



To include the radiative cooling processes from the heavy elements, we use the method 
described in Smith, Sigurdsson, & Abel (2007), in preparation. The nonequilibrium abun- 
dances and cool ing rates of H, H "*", H~, He, He"*" , He"*"^ , H2, H^, and e~ are calculated 
internally, as in I Abel et al.l (120021 ): lAnninos et al.l (119971 ). Meanwhile, the metal cooling 
rates are interpol ated from large grids of values, precomputed with the photoionization soft- 
ware, CLOUDY (IFerland et al.lll998l ). We ignore the cooling from dust and focus only on 
the contribution of gas-phase metals in the optically-thin limit. Unlike other studies of the 
formation of the first metal-enriched structures, we do not assume the presence of an ionizing 
UV background. In our model, the singular pop III star that was associated with the dark 
matter halo in which our stars form has already died in a supernova. We also assume any 
other Pop III stars are too distant to affect the local star-forming region and that QSOs have 
yet to form. We use the coronoal equilibrium command when constructing the cooling 
data in CLOUDY to simulate a gas where all ionization is coUisional. The metal cooling 
data was created using the Linux cluster, Lion-xo, run by the High Performance Computing 
Group at The Pennsylvania State University. As a consequence of our choice to ignore any 
external radiati on, we do not observe th e fine-structure emission of [C ll] (157.74 fim) that 
was reported by ISantoro fc Shull J2006h to be important. Instead, cooling from C comes in 
the form of fine-structure lines of [C l] (369.7 /im, 609.2 /^m). The cooling from [C l] in our 
st udy dominates iri the sa me range of densities and temperatures as the cooling from [C ll] 

We observe the contributions of the other coolants studied by 



m 



Santoro fc Shu. 



3 J2006h 



Santoro &: Shulll (120061 ) . [O l], [Si ll], and [Fe ll], to be in agreement with their work. In 



addition, we find that emission from [S l] (25.19 /im) dominates the cooling from metals at 
n ~ 10^ cm~^ and T ~ 1-3 x 10^ K. The absence of UV radiation in our simulations also 



allows H2 to form, differentiating this study from iBromm et al.l (120011 ). This allows for a 
more direct comparison between the metal-free and metal-enriched cases. 

The simulations are run until one or more dense cores form at the center of the dark 
matter halo and a maximum refinement level of 28 is reached for the first time, giving us 
a dynamic range of greater than 10^°. Only the simulation with Z = 10~^ Zq reached 28 
levels of refinement. The three other simulations were stopped after reaching 27 refinement 
levels, since their central densities were already higher than the simulation with Z = 10~^ 
Zq. Table 1 summarizes the final state of each simulation, where Zcoi is the collapse redshift, 
Imax is the highest level of refinement, rimax is the maximum gas density within the box, and 
Atcoi is the time difference to collapse from the metal-free simulation. 
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3. Results 

As can be seen in Table 1, the runs with higher metallicities reach the runaway collapse 
phase faster. The relationship between metallicity relative to solar and Atcoi is well fit by 
a power-law with index, n ~ 0.22. Gas-clouds with more metals are able to radiate away 
their thermal energy more quickly, and thus, collapse faster. An inverse relation between 
metallicity and the number of grids and grid-cells exists because the low-density, background 
gas evolves at roughly the same rate in all simulations, yet has more time, in the runs with 
lower metallicities, with which to collapse to higher density, requiring additional refinement. 

Our simulations, shown in Figure 1, display a qualitative transition in behavior between 
metallicities of 10~^ Zq and 10"'^ Zq. In the runs with the highest metallicities (Figure IC 
and ID), the central core is extremely asymmetric, and multiple density maxima are clearly 
visible. All four runs display similar large-scale density profiles (Figure 2A). Radiative cooling 
from H2 becomes extremely inefficient below T ~ 200 K, creating the effective tempe rature 
floor, visible in Figure 2B for the metal- free case (lAbel et al.ll2002l : iBromm et al.ll2002l ). At n 
~ 10^ cm~^, the rotational levels of H2 are populat ed according to LTE, reducing the coo ling 
efficiency and causing the temper atur e to in crease (lAbel et al.ll2002l : iBromm et al.ll2002l ). In 
the isothermal collapse model of IShul (Il977l ). the accretion rate is proportional to the cube 
of the sound speed. The increase in temperature leads to an increase in the accretion rate, 
causing the density, and thus, the enclosed mass (Figure 2C), to be slightly higher inside the 
central ~ 0.1 pc in the metal-free case. A similar situation occurs further within for the Z 
= 10~^ Zq and, later, the 10~^ Zq cases, as the metal cooling is overwhelmed by adiabatic 
compression heating and the temperature begins to rise with density. The presence of metals 
at the level of 10~^ Zq enhances the cooling enough to lower the gas temperature to ~ 75 
K. Metallicities greater than 10~^ Zq provide sufficient cooling to bring the gas down to 
the temperature of the cosmic microwave background, where Tcmb ^ 2.7 K (1 + z ) . The 
gas temperatures are in general agreement with the calculations of lOmukai et al.l (120051 ) 
that include a CMB spectrum at z = 20. Fragmentation requires that the cooling time be 
less than the dynamical time. Figure 2D shows that this criterion is essentially never met 
in the zero metallicity case, and only marginally in the Z = 10"'^ Zq case. However, the 
fragmentation criterion is more than satisfied in the Z = 10~^ Zq and 10~^ Zq cases over a 
wide mass-range. 



on 



In order to locate f ragments within our simulations, we employ an algorithm, based 



Williams et al.l (119941 ). that works by identifying isolated density countours. Before we 



search for clumps, we smooth the density field by assigning each grid-cell the mass-weighted 
mean density of the group of cells including itself and its neighbors within one cell-width. 
This serves to eliminate small density perturbations that would be misidentified as clumps by 
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the code. In order to directly compare the fragmentation from each simulation, we limit the 
search for clumps to the 1 Mq of gas surrounding the cell with the highest density. On larger 
scales, all of the runs display a filamentary structure that is qualitatively similar. No other 
region in any of the simulation boxes has collapsed to densities comparable to those found 
within the region where the clump search is performed. The results are shown in Figure 3. 
A single clump exists in the metal-free and 10~^ Zq simulations, containing 99.7% of the 
total mass within the region of interest. In the simulation with Z = 10~^ Zq, 91% of the 
mass is shared between two clumps with 0.52 Mq and 0.39 Mq. In the same simulation, we 
also find two smaller clumps 0.06 Mq and 0.02 Mq. Finally, in the Z = 10~^ Zq simulation, 
we see two clumps with 0.79 Mq and 0.21 Mq. 



4. Discussion 



We have shown, through three-dimensional hydrodynamic simulations, that fragmen- 
tation occurs in collapsing gas with metallicities, Z > 10^'^ Zq. Our results indicate that 
star- formation occurs in exactly the same manner at metallicity, Z = 10^^ Zq, as it does at 
zero metallicity. The similarities between the simulations with metallicities, Z = 10~^ Zq 
and 10~^ Zq, suggest that the transition to low-mass star-formation is complete by 10"^ Zq, 
implying that the entire transition occurs over only one order of magnitude in metal abun- 
dance. More simulations, bracketing the metallicity range, 10"^ to 10-3 Zq, will test how 
abrupt the transition truly is. We will also explore the effect of non-solar abundances on the 
low metallicity IMF. It has been recently argued that dust cooling at high densities (n > 10^ ^ 



cm 3) can induce fragm entation for rnetalli cities as low as 10 ^ Zq ISchneider et al.l (120061 ). 



In light of the work by iFrebel et al.l (120071 ). who note the absence of stars with Dfrans < 
-3.5, where Dtrans is a measure of the combined logarithmic abundance of C and O, it seems 
unlike ly that Zr. r is this low. While the fragmentation mode discussed in ISchneider et al. 



( 120061 ) ■ and also lOmukai et al.l (120051 ) . may truly exist, it is possible that metal yields from 
Pop III supernovae overshoot this metallicity, for realistic mixing scenarios, leaving almost 
no star -forming regions with such a low concentration of heavy elements. Similar to our 
results, lOmukai et al.l (120051 ) note that only high-mass fragments are produced when Z = 
10"^ Zq. If Pop III supernovae are able to immediately enrich the local universe to Z = 10"^ 
Zq, the high-density dust cooling fragmentation mode would be skipped altogether, and the 
high-mass stars that formed via the mode observed at 10^^ Zq would leave no record in the 
search for low-metallicity stars in the local universe. 



We have limited the search for fragments to the dense 1 Mq core at the center of each 
simulation. Within this region, it is unlikely that any more fragments will form in any of 
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the simulations. In all of the cases presented, the cooling has begun to be overwhelmed 
by compression heating such th at the ce i itral t emperature is now increasing with increasing 
density, which was indicated by iLarsoru (120051 ) to be the end of hierarchical fragmentation. 
Fragmentation may continue in the surrounding lower density gas in the cases of Z = 10"'^ 
Z0 and 10~^ Zq. The final stellar masses of these objects will also be affected interaction 
and accretion that will occur in later stages of evolution. In the two lowest metallicity cases, 
the gas immediately surrounding the central core evolves slowly enough that it will not have 
sufficient time to reach high densities befo re the UV radi a tion f rom the central, massive star 
dissociates all of the H2. As was shown by lBromm et al.l (I2OOII ). clouds with metallicities, Z 
< 10~^ Z0 are unable to collapse without the aid of H2 cooling. 

In the two simulations in which significant fragmentation is obse rved, Z = 10~^ and 



10 ^ Z0, the gas is able to cool rapidly to the temperature of the CMB. IWise fc Abel! (120051 ) 



predict that the rate of Pop III supernovae peaks at a redshift, z ~ 20, and then drops 
off sharply, implying that metal production from Pop III stars is effectively finished at this 
point. In this epoch, the characteristic r nass-scale for metal-en riched star formation will be 
regulated by the CMB, as is predicted in lBromm fc Loebl (120031 ). Thus, the first Pop II stars 
wi ll be considera bly more massive, on average, than stars observed today, as was suggested 
by iLarsonl (119981 ). Observations of low - mass preste llar cores in the local universe reveal them 
to have temperatures of about 8.5 K ( Evanslll999l ). implying that the IMF may not become 
completely 'normal' until z < 3 when the CMB fell below this temperature. 
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Table 1 

Simulation Final States 

Z (Zq) Zcoi Kax Grids Cells Umax (cm~^) Atco; (Myr) 

18.231519 27 8469 4.82 xlO^ 4.11 xlO^^ 

10-^ 18.838816 27 8060 4.64 xlO^ 3.90 xlO^^ 9.19 

10-=^ 19.336557 27 7911 4.56 xlO^ 1.65 xlO^^ 16.21 

10-2 20.032518 28 7521 4.42 xlO^ 1.50 xlO^^ 25.33 



Fig. 1. — Slices through gas density for the final output of simulations with Z = (A), 10~^ 
Zq (B), 10~^ Z0 (C), and 10~^ Z0 (D). Each shce intersects the grid-cell with the highest 
gas density and has a width of 2 x 10~^ of the computation box, corresponding to a proper 
size of ~4 X 10~^ pc (84 AU). The color-bar at bottom ascends logarithmically, from left to 
right, spanning exactly four orders of magnitude in density. 
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radius [AU] enclosed gas mass [M^] 




radius [pc] enclosed gas mass [MJ 



Fig. 2. — Radially averaged, mass-weighted quantities for the final output each simulation: 
Z = (red), 10~^ (green), 10~^ Z© (blue), and 10~^ Z0 (purple). A: Number density vs. 
radius. B: Temperature vs. enclosed mass. C: Enclosed gas mass vs. radius. D: Ratio of 
crossing time to cooling time vs. enclosed mass. The classical criterion for fragmentation is 
met when the ratio of the crossing time to the cooling time is greater than 1. 
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Fig. 3. — Masses of clumps found within the final output of each simulation. The location 
on the X and y axes corresponds to the log of the clump mass and the metallicity of the 
simulation. Colors are the same as in Figure 2. The radii of the circles are proportional to 
the masses of the clumps they represent. A factor of 10 in mass is equivalent to a factor of 
2 in radius. The search for clumps is limited to the 1 Mq surrounding the grid cell with the 
highest gas density. Only clumps with at least 1000 cells are plotted. 



